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Abstract 


A commonly  encountered  integer  linear  program,  basic 
to  cyclic  staffing  and  scheduling,  has  a constraint  matrix 
possessing  the  property  of  "circular  l's  in  columns."  In 
general,  such  a matrix  is  not  unimodular,  balanced,  or  per- 
fect. Nevertheless,  many  such  problems  may  be  efficiently 
solved  for  integer  answers.  A change  of  variable  transforms 
them  to  comfortably  finite  and  reassuringly  predictable 
series  of  minimum  cost  network  flow  problems. 
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CIRCULAR  l's  AND  CYCLIC  STAFFING 


1.  Two  Fundamental  Staffing  Models 

Consider  the  integer  linear  program 

min  cx 

s . t . Ax  > b (1.1) 

x > 0,  integer 

where,  throughout  the  paper,  b and  c are  vectors  with  all 
entries  integer  and  A is  an  m x n matrix  with  all  entries 
0 or  1.  Without  loss  of  generality,  we  may  assume  b,  c > 0. 

To  represent  continuous  workshifts  in  linear  time,  a 
common  staffing  model  has  A possess  the  property  of  conse- 
cutive l's  in  columns  (e.g.,  Veinott  and  Wagner  [19]). 

Such  matrices  are  happily  met  since  they  are  known  to  be 
totally  unimodular;  moreover,  for  such  matrices,  problem 
(1.1)  is  equivalent  under  linear  transformation  to  the 
minimum-cost  network  flow  problem 

min  cx 

s.t.  [TA,  -T] x = Tb  (1.2) 

x > 6,  integer 

where  T is  the  m x m matrix 


1 


(1.3) 


and  where  "equivalence"  means  here  that  x solves  (1.1)  if 
and  only  if  x solves  (12)  [11,12].  Transforming  (l.i)  by 

T to  reveal  its  network  structure  corresponds  to  succes- 
sively row-reducing  the  constraints  of  (1.1)  [ 19 ] . 

Since  the  minimum  cost  network  flow  algorithm  is  formally 
efficient  [8  ],  we  may  consider  (1.1)  to  be  efficiently 
solvable  in  its  guise  (1.2). 

The  second  basic  staffing  model  represents  continuous 
workshifts  in  cyclical  time  [3].  For  this  model  the  matrix 
A possesses  the  property  of  circular  l's  in  columns  [18], 
as  for  instance  in  Example  1.1,  where  the  strings  of  l's 
may  be  imagined  to  wrap  around  the  matrix.  Such  matrices 
are  in  general  neither  unimodular,  balanced,  nor  perfect  [16]. 
Indeed  they  are  notorious  for  the  fractional  extreme  points 
which  they  induce  in  (1.1)  [13]. 

A special  n x n circular  l's  matrix  has  in  each  column 
a band  of  k l's  permuted  cyclically  (see  Examples  1.2  and 
1.3).  We  will  call  these  (k,  n)  matrices. 

The  most  fundamental  of  the  cyclic  staffing  models  is 
given  by 

min  lx 

s . t . Ax  > b (1.4) 

x > 0,  integer 
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110  11 
110  0 0 
0 10  0 0 
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Example  1.1: 


A matrix  with  "circular  1 


s in  columns. " 


110 
10  1 
Oil 


10  0 11 
110  0 1 
1110  0 
0 1110 
0 0 111 


Example  1. 2 : A (2,  3)  matrix.  Example  1.3:  A (3,  5)  matrix. 
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where  A is  a (k,  n)  matrix.  The  objective  corresponds  to 
minimizing  the  total  workforce  size  necessary  to  meet 
period  manpower  requirements  b.  For  A a (5,  7)  matrix, 
this  problem  was  studied  by  Tibrewala,  Philippe,  and  Browne 
[17]  (also  by  various  others  [2,4,  6,15]),  where  the  (5,  7) 
matrix  represents  workshifts  with  two  consecutive  days  off 
each  week.  They  observe  that  their  solution  generalizes 
to  A a (k,  k + 2)  constraint  matrix.  A rather  more  complex 
solution  technique  is  proposed  by  Guha  [10]  for  general  (k,  n) 
matrices.  In  this  study  we  generalize  problem  (1.4)  in  two 
ways  and  offer  considerably  simpler  and  formally  efficient 
solutions . 


2 . Transformations  of  Variables 
Consider  the  problem 


min  cx 

s . t . Ax  > b (2.1) 

x > 0 


which  we  may  write  as 


min 

cx 

s . t . 

‘ A ‘ 

" b • 

X > 

ft 

L 1 J 

0 

(2.2) 


Now  let  T be  a nonsingular  matrix  and  consider  the  change  of 
variables  x = Ty.  Since  T is  nonsingular,  (2.2)  is  equiv- 
alent to 
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min  (cT)y 


s.  t. 

1 

> 

1 

fbl 

L T_ 

A 1 

>1 

i 

IO 

1 

y unrestricted 

in  the  sense  that  if  x is  feasible  to  (2.2),  then  T is 
feasible  to  (2.3),  and  if  y is  feasible  to  (2.3),  then  Ty 
is  feasible  to  (2.2).  If  in  addition  T is  unimodular,  then 
if  x has  all  integer  entries,  T 1x  has  all  integer  entries, 
and  if  y has  all  integer  entries  Ty  has  all  integer  entries. 
Therefore , 

Observation  2.1;  For  T nonsingular  and  unimodular  the  in- 
teger-constrained versions  of  problems  (2.2)  and  (2.3)  are 
equivalent  in  the  sense  that  x solves  (2.2)  iff  y = T ^x 
solves  (2.3) . 

We  will  use  this  insight  to  construct  equivalent  integer 
programs  wherein  special,  exploitable  structure  is  displayed. 

3 . Almost  a Network 

A key  idea  of  this  paper  is  that  under  certain  condi- 
tions, when  A is  a circular  l's  matrix,  problem  (1.1)  is 
"almost"  a network  flow  problem.  Problem  (3.1)  below, 
where  A is  a (3,  5)  matrix,  will  provide  a continuing  illus- 
tration of  this  class  of  problems.  Later  we  will  observe 
that  the  ideas  to  follow  generalize  easily. 


r 


Moreover,  since  such  a T is  both  nonsingular  and  unimodular, 
(3.2)  is  equivalent  to  (3.1)  as  an  integer  linear  program. 
Thus  solving  (3.2)  solves  (3.1).  We  solve  (3.2),  on  the 
strength  of  the  following, 

Observation  3.1:  With  the  exception  of  the  last  column, 
that  corresponding  to  y,.,  problem  (3.2)  is  the  linear  pro- 
gramming dual  of  a network  flow  problem. 

That  is,  if  we  fix  (temporarily)  y^,  problem  (3.2)  may 
be  written  as 


yl'  y2'  y3'  y4  unrestricted  but  integer 


which  is  the  linear  programming  dual  of  a network  flow 
problem.  The  obvious  idea  is  to  fix  values  of  y,.  over  its 
allowable  range  and  solve  corresponding  network  flow  pro- 
blems until  the  best  objective  value  is  found.  Let  us  re- 
fine and  extend  this  idea. 

7 
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4 . Properly  Compatible  l's 

Following  Tucker  [18] , we  define  a 0 - 1 matrix  A 

to  have  properly  compatible  circular  l's  in  columns  if 

and  only  if  (i)  the  l's  in  each  column  are  circular,  and 

(ii)  for  any  two  columns  a.  and  a,  , if  the  first  (in  a 

J K 

cyclic  sense)  1 in  a . preceeds  that  of  a.  , then  the 

D K 

last  (in  a cyclic  sense)  1 in  ak  does  not  preceed 
that  of  a...  Roughly  speaking,  if  a circular  band  starts 
later  than  another,  it  can  end  no  earlier.  The  matrices  of 
Example  4.1  illustrate  properly  compatible  circular  l's  in 
columns.  For  matrices  with  this  property,  a natural 
ordering  of  the  columns  suggests  itself. 

Column  Orderinq  Alaorithm 


1.  Order  columns  in  groups,  where  group  i consists 
of  those  aj  whose  first  (in  a cyclic  sense)  1 appears  in 
row  i.  Then, 

2.  Within  each  group,  order  columns  so  that  aL  pre- 
ceeds a,  if  the  last  (in  a cyclic  sense)  1 of  a . preceeds 

K D 

that  of  a^.  The  columns  of  Example  4.1  have  been  so 
ordered.  Henceforth,  we  assume,  without  loss  of  generality, 
that  a matrix  with  properly  compatible  l's  in  columns  has 
its  columns  ordered  as  above.  Important  for  us  shortly  will 


Observation  4.1:  A matrix  with  properly  compatible  circular 
l's  in  columns  has  the  property  of  circular  l's  in  rows. 


Consider  now  the  problem 
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Properly  compatible  circular  l's  in  columns. 


Example  4,lb;  A circular  l's  matrix,  the  columns  of  which 
are  not  properly  compatible. 
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min  cx 

s . t . Ax  > b (4.1) 

x > 0 , integer 

where  A has  properly  compatible  l's  in  columns.  Perform 
the  change  of  variables  given  by  x = Ty,  where  T is  the  non- 
singular unimodular  matrix  defined  in  (1.3).  Then,  we  have 
an  equivalent  integer  linear  program 


min  (cT)y 


;.t.  r at  i r b " 

L T J y - L 5 . 


(4.2) 


y unrestricted  but  integer 


Since  A has  circular  l's  in  rows,  each  row  r^  of  A has 
either  consecutive  l's  or  consecutive  0's  [18].  Therefore 


is  of 

the 

form 

(i) 

?i 

= (0 

(ii) 

r . 

= (1 

l 

(iii) 

?i 

= (1 

(iv) 

?i 

= (0 

But  then  each  row  r^T  of  AT  is  of  the  form 

(i)  r^T  ~ (0,...,0,  — 1,  0,...,0,  1,  0,*..,0),  or 

(ii)  r^T  = (0,...,0,  1,  0 , • • • , 0 , —1,  0,.«»,0,  1),  or 

(iii)  r^T  = (0,...,0,  1,  0,...,0),  or 

(iv)  r /r  = (0 , . . . , 0 , -1,  0 , . . . , 0 , 1),  respectively. 


Note  that,  excluding  the  n*"*1  column,  each  row  of 


AT 

T 


has  at 
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most  one  +1  and  one  -1,  all  other  entries  being  0. 

For  notational  convenience,  let  us  partition  T into 
its  nth  column  and  the  remainder  of  the  matrix:  T * [Tr,  tnl 
= [Tr,  enl , since  the  nth  column  of  T is  * (0,...,0,  1). 
Similarly  partition  y = (yr»  yfi)  . Then  problem  (4.2)  may 
be  rewritten  as 


min  (cTr)yr  + cnyn 


" at 

• “ 

a 

y 

’ b “ 

r 

n 

r 

> 

_ Tr 

®n  _ 

_ yn  J 

0 

yr,  yn  unrestricted  but  integer 


Now  we  can  formally  state 

Lemma  4.1:  If  for  problem  (4.1),  A has  properly  compatible 
circular  l's  in  columns,  then  under  the  prescribed  change  of 


variables 


AT, 


is  the  transpose  of  a network  matrix. 


That  is,  for  fixed  integral  yn,  the  resultant  version  of  (4.3) 


min  (cTr)yr 


u 

b* 

< 

i 

’ E - v. ' 

_ Tr_ 

yr  - 

. 5 ' Vn  . 

yr  unrestricted  but  integer 


(4.4) 


is  the  linear  programming  dual  of  a network  flow  problem. 

Thus,  problem  (4.4)  is  efficiently  solvable,  at  least  through 

its  dual.  This  suggests  the  idea  of  searching  through  the 

allowable  values  of  yn,  solving  a tractable  subproblem  (4.4) 

each  time,  to  find  a (y*,  y*)  which  minimizes  (cT  )y  + c y . 

r n r r n n 


ii 


5.  Stalking  the  Wild 

First  we  determine  the  allowable  range  of  the  integer 

yn.  Let  y*  be  a value  of  yn  in  some  optimal  solution  to 

(4.2)  » and  let  bmax  be  the  largest  entry  in  b.  Then 

Lemma  5.1:  b < y*  < lb  for  some  y* 

— max  - n - n 

Proof : Since  y = T-1x,  and 

1 

1...1  . 


we  have  that  yfi  = lx.  To  show  the  lower  bound,  it  is 


sufficient  to  observe  that  sihce  x > 0,  y = lx  > Ea. .x.  > 

- -'n  - j i]  ] - 

b.  Vi.  Therefore,  y*  > 

To  establish  the  upper  bound,  we  may  assume  that  at 
optimality  every  variable  in  (4.1)  appears  in  some  tight 
constraint,  since  otherwise  that  variable  could  be  reduced, 
feasibility  maintained,  and  the  objective  function  not  in- 
creased. Summing  over  the  set  S of  tight  constraints  yields 
E 


Ed.  .x*  = E b.  ; but  lb  _>  E b . = E E d .x!  > = y 

if S j 13  3 itS  1 ieS  1 ieS  j 3 3 


Q.E.D. 

To  remind  us  of  its  dependence  on  yn,  let  the  objective 

function  of  problem  (4.4)  be  written  as  (cTr)yr  = z (yn) 

and  let  the  optimal  value,  for  fixed  yn,  be  z*  (yn) . 

Lemma  5.2:  z*(y_)  is  convex  in  y„  over  b < y_  < 15. 

" n it  iuax  * n " 

Proof:  Since  the  constraint  matrix  of  problem  (4.4)  is 

totally  unimoduLar,  the  integral  restrictions  may  be  dropped. 
Then  the  desired  conclusion  follows  from  similar  results 
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for  continuous-valued  linear  programs  (e.g.,  Geoffrion  and 
Nauss  [ 9 ] ) . 


Q.E.D. 

Lemma  5.3:  The  optimal  function  value  (cT)y  of  problem 
(4.2)  is  convex  in  y over  b < y < lb. 

Proof : Clearly  cnyn  is  convex  in  yn;  z* (yn)  is  convex  in 

1'.  by  Lemma  5.2,  so  since  sums  of  convex  functions  are  convex 
i ) + c^y  . is  convex  in  y^.  But  this  is  the  optimal 
function  value  of  problem  (4.3)  and  therefore  of  problem 
(4.2) 

Q.E.D. 

6.  A Solution  Technique 
Given  the  problem 


min  cx 
s . t . Ax  b 

- = . . (6.1) 

x ^ 0,  integer 

where  A has  properly  compatible  circular  l's  in  columns, 
the  preceding  results  justify  this  solution  procedure, 

Step  0:  Perform  the  change  of  variables 

Let  x = Ty,  where  T is  defined  by  (1.3)  to  form  the  equivalent 
problem 


min  (cT)y 


— — 

— — - 

s . t . 

AT 

b 

y > 

T 

0_ 

(6.2) 

y unrestricted  but  integer 

Step  l:  Solve  the  equivalent  problem , (6.2) 


A.  Note  bounds  on  integer  y*:  b < v*  < lb. 

B.  Minimize  z*  (yn ) + cnyn  = (cT)y  over  this  interval. 

Since  y is  integer  and  z* (y  ) + c y is  convex  in  y , 
n n n-1  n 1 n 

an  efficient  technique  such  as  Fibonacci  search  [21]  may  be 
used.  Furthermore,  for  fixed  y^,  z* (yn)  is  readily  calcul- 
ated by  solving 

z* (yn)  = min  (cTr)yr 


b - a y 
n n 

yr  - 0 - e y 

nJn 


yf  unrestricted  but  integer 


(6.3) 


Since  this  is  the  dual  of  a minimum  cost  network  flow  prob- 
lem, it  is  efficiently  solvable.  Let  (cT)y*  = min  z*  (y^)  + c^y^ 

yn 

and  let  (y*,  y*)  = y*  be  the  associated  solution;  then  y* 
solves  (6.2)  and  (cT)y*  is  the  optimal  function  value. 

Step  2:  Construct  the  optimal  solution  to  (6.1)  by 
the  change  of  variables  x*  = Ty*. 


7 . Efficiency  of  the  Algorithm 

This  solution  procedure  works  efficiently,  even  for 
pessimists,  by  the  following  worst-case  analysis. 

Step  0,  the  initial  change  of  variables,  requires  no 
more  than  0(mn)  steps. 

Step  1 requires  the  solution  of  (6.3)  for  fixed  y^. 

But  the  network  flow  algorithm  solves  (6.3)  in  a number  of 
steps  which  is  bounded  above  by  a polynomial  in  the  size  of 
the  encoding  of  the  problem  data  [8].  We  may  take  this 


polynomial  to  be  p(m,  n,  log2  lb,  log2  lc,  Yn) • But  since 
y*  < lb,  log-y*  <_  log-  lb,  so  that  we  may  consider  the  solu- 
tion  to  (6.3)  to  require  no  more  than  0(p(m,  n,  log2  lb, 
log2  lc) ) for  some  polynomial  p.  And  since  Fibonacci  search 
requires  that  we  consider  no  more  than  0 (log^  lb)  values  of 
y , we  may  determine  y*  = (y*,  y*)  is  no  more  than  Odog^lb 
•£(m,  n,  log2  lb,  log2  lc))  steps. 


Step  2,  the  final  change  of  variables,  requires  0(n) 


steps . 


Therefore  the  solution  procedure  solves  (6.1)  in  at 
worst  0 (mn  + log3  lb  • p(m,  n,  log2  lb,  log2  lc) ) steps. 
Since  this  is  polynomial  in  a binary  encoding  of  the  pro- 
blem data  [1],  we  have  proven 

Lemma  7.1;  Problem  (6.1)  is  solved  by  the  solution  tech- 
nique with  formal  efficiency  relative  to  a binary  encoding 
of  the  problem  data. 


8 . A Special  Objective  Function 

For  a special  objective  function,  a wider  class  of 
problems  may  be  solved  and  additional  results  discovered. 


Consider 


min  lx 


s . t . Ax  > b 


(8.1) 


x > 0,  integer 


where  A displays  the  property  of  circular  l's  in  columns 
(not  necessarily  properly  compatible) . 
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x*  = (x*, . . . ,x* , . . . ,x*, . . . ,x*)  solve  (8.1).  Then 

a j k n 

(x*,...,x*  + x*,  . . . ,0.  , . . . ,x*)  is  feasible  to  (8.1)  and, 

-L  3 ^ ^ n 

moreover,  has  the  same  (optimal)  objective  value.  Therefore 

Lemma  8.1:  An  optimal  solution  to  problem  (8.1)  exists 

for  which  none  of  the  columns  of  A corresponding  to  nonzero 

variables  are  dominated  by  any  other  such  column  of  A. 
Therefore  we  may  reduce  (8.1)  by  eliminating  any  columns 


of  A (and  associated  variables)  which  are  dominated.  But 
then  the  resulting  matrix  displays  properly  compatible  cir- 
cular l's  in  columns,  so  that  the  problem  is  solvable  by 
the  approach  just  presented.  (Note  that,  in  fact,  it  is 
sufficient  for  this  conclusion  to  assume  so-called  "agree- 
able" costs,  for  which  c . < c,  iff  a . > a,  (cf.,  [14])). 

1 - K 1 ~ k 

Let  us  assume  that  the  matrix  A has  been  pruned  of 

dominated  columns.  Then  the  special  properties  of  the 

transformed  problem  are  of  interest.  In  particular,  the 

new  objective  function  is  (cT)y  = eny  = 0yr  + yn-  Thus 

solving  the  transformed  problem  (6.2)  is  tantamount  to 

finding  the  smallest  integer  yn  for  which  the  constraints 

of  (6.3)  have  a feasible  solution.  Equivalently,  we  seek 

the  smallest  integer  y for  which  the  dual  network  flow 

n 

problem  to  (6.3)  is  not  unbounded,  i.e.,  is  free  of  cycles 
of  positive  net  cost. 

For  the  special  objective  function  cx  such  that 
C1  - c2  - *•*  - Cn-1  includes  the  objective  function 

lx) , a particularly  simple  solution  technique  applies  to  the 
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transformed  problem  (6.3).  The  new  objective  function  has 
the  property  (cT^)  > 0;  furthermore 

one  +1  in  each  row,  and  at  least  one  +1  in  each  column.  Thus 
this  version  of  (6.3)  is  solvable  by  the  simple  recursive 
substitution  scheme  of  Dorsey,  Hodgson,  and  Ratliff  [ 7 ]. 


has  no  more  than 


9.  Close  Enough 

For  the  transformed  version  of  problem  (8.1),  an 
interesting  round-off  result  holds  (see  similar  results  in 
[4,20]).  Recall  that  the  transformed,  equivalent  version  of 
(8.1)  is 


min  Oy  + y 
. 1 r Jrn 


” AT 

a 

y 

" b ' 

r 

n 

Jr 

> 

T 

e 

y 

0 

r . 

n 

L n _ 

yr • yn  unrestricted  but  integer 


Lemma  9.1:  Let  y'  = (y^,...,y^)  solve  the  continuous-valued 
relaxation  of  (9.1)  then  y*  = ( Tyil  * * - - - * > solves 

the  integer-restricted  problem  (9.1). 

Proof:  Clearly  r»;i  is  a lower  bound  on  the  optimal  function 

value  of  (9.1).  Moreover  ( fy^l  • [V^l  • ' ' ' • fy^l  ) is  an  integer- 
valued vector  which  achieves  this  value.  To  see  that  this 
vector  is  feasible  to  (9.1),  we  will  show  that  it  satisfies 
each  of  the  constraints,  of  which  there  are  three  types: 


(i) 

yj 

- *k 

i bi 

(ii) 

yj 

- yk 

+ *n  i 

(iii) 

yj 

•H 

A 

A 1 
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First  observe  that  for  any  two  numbers  a and  b. 


- Tal  = f-aj  , (9.2) 

and 

fa"j  + fb"|  > Ta+t)l  (9.3) 

By  (9.3)  , f a-b]  + fb"J  > fa"]  so  that  M - M > -Ml  - Lb-aJ 
by  (9.2).  Then  by  the  last  inequality  we  have 

(i)  l~y jl  - [v£\  ^ |_yj-ykJ  - LbiJ  = bi  since  bi 

is  integer. 

(ii)  [y]\  - ry’l  + Ty^l  > Ty^+y;!  - [y^\  > 

|yj-y^+y^J  > ^ij  = bi  since  bi  integer. 

dii)  f y jl  1 Y]  1 v 

Hence,  each  of  the  constraints  of  (9.1)  is  satisfied  and 
(fy-],^'! an  opti"131  feasible  solution. 

Q.E .D. 

Therefore  problem  (8.1)  may  be  solved  by  the  following 
simple  application  of  linear  programming: 

(i)  Solve  the  continuous-valued  relaxation  of 
(8.1)  by,  for  example,  the  simplex  method  of  linear  program- 
ming. Let  the  solution  be  x’ . 

(ii)  Transform  the  solution  via  y'  = T ^x' , for  T 
as  in  (1.3). 

( iii)  Round-up  y*  = ( fy}"|  , fy^]  * • • • * [y^\  > • 

(iv)  Transform  back  to  x*  = Ty*.  Then  x*  solves 
the  integer  program  (8.1). 
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10. 


10.  Applications 

A.  Cyclic  Staffing  with  Over*  time 

A basic  staffing  problem  involves  a facility  such  as  a 
hospital  that  operates  24  hours  each  day.  Assume  there  are 
fixed  hourly  staff  requirements  b^ , and  that  there  are  three 
basic  work  shifts,  each  of  eight  hours  duration:  0700-1500, 
1500-2300,  and  2300-0700.  Overtime  of  up  to  an  additional 
eight  hours  is  possible  for  each  shift.  What  is  the  minimum 
cost  number  of  workers  and  their  shifts  such  that  all  staff 


requirements  are  met?  This  problem  may  be  formulated  as 
in  Figure  10.1,  where  the  constraint  matrix  displays  pro- 
perly compatible  circular  l's  in  columns.  Thus  the  problem 


is  efficiently  solvable  by  a bounded  series  of  network  flow 
problem. 

B.  Days-off  Scheduling 

A problem  studied  by  Brownell  and  Lowerre  [5]  is  to 


minimize  the  total  workforce  necessary  to  meet  daily  staff- 
ing requirements,  where  each  worker  is  guaranteed  two  days 
off  each  week,  including  every  other  weekend.  For  the  case 
in  which  the  days  off  each  week  are  to  be  consecutive,  the 
problem  may  be  formulated  as  in  Figure  10.2.  The  rows  of 
the  matrix  display  more  complicated  cyclic  structure  than 
simple  circular  l's;  but  since  the  matrix  has  circular  l's 
in  rows,  the  same  change  of  variables  transforms  the  problem 
to  efficiently  solvable  form. 


min  cx 


s.t. 

07  [ 1 1 1 1 1 1 1 1 1 000000000  Ollllllll" 

08  111111111  000000000  001111111 

09  111111111  000000000  000111111 

10  111111111  000000000  000011111 

11  111111111  000000000  000001111 

12  111111111  000000000  00C000111 

13  111111111  000000000  000000011 

14  111111111  000000000  000000001 

15  011111111  111111111  000000000 

16  001111111  111111111  000000000 

17  000111111  111111111  000000000  - r 

18  000011111  111111111  000000000  x- 

19  000001111  111111111  000000000 

20  000000111  111111111  000000000 

21  000000011  111111111  000000000 

22  000000001  111111111  000000000 

23  000000000  Ollllllll  111111111 

24  000000000  001111111  111111111 

01  000000000  000111111  111111111 

02  000000000  000011111  111111111 

03  000000000  000001111  111111111 

04  000000000  000000111  111111111 

05  000000000  000000011  111111111 

06  [ 0 0 0 0 0 0 0 0 0 000000001  111111111 


x > 0,  integer 


Figure  10.1:  A cyclic  staffing  problem  with  overtime. 
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